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Abstract 

The Polyakov world-line path integral describing the propagation of gluon field 
quanta is constructed by employing the background gauge fixing method and is 
subsequently applied to analytically compute the divergent terms of the one (glu- 
onic) loop effective action to fourth order in perturbation theory. The merits of the 
proposed approach is that, to a given order, it reduces to performing two integra- 
tions, one over a set of Grassmann and one over a set of Feynman-type parameters 
through which one manages to accommodate all Feynman diagrams entering the 
computation at once. 
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1 Introduction 



Improved methods, in comparison to the Feynman diagrammatic ones, for expediting 
perturbative calculations in QCD have emerged, within the last decade or so, through 
the adoption of first-quantization-based approaches. The latter involve either string or 
world-line agents through which one describes the field theoretical system. The original 
efforts in this direction were string inspired and were based on realizations, made in 
the late 1980's [1-5], regarding the relation between string and non-abelian gauge field 
theories in the limit of an infinite string tension. Following their own involvement in 
such studies, Bern and Kosower [6] established a set of rules expediting efficient, one 
loop, computations in non-abelian gauge theories. Through them one could encompass 
contributions of a host of Feynman diagrams at once. Further extensions of the string- 
inspired approach were subsequently carried out in Refs. [7-12]. 

World-line based methodologies aiming at the same goals soon followed through the 
work of Strassler [13] who proposed suitably defined (one-dimensional) path-integrals 
for the various quantum field systems he considered. Extensive use was made of su- 
persymmetric one-dimensional particle coordinates [14-16] 'living' on the paths. The 
end result was the accomplishment of the full reproduction of the Bern-Kosower rules. 
Strassler's approach was further pursued in Refs. [17-19], where computations pertain- 
ing to multi-loop configurations in QED as well as effective actions involving constant, 
external (chromo)electric and (chromo)magnetic fields were undertaken. 

Now, the world-line casting of relativistic quantum systems is an old story, which 
goes back to Fock [20], Feynman [21] and Schwinger [22]. Notable, relevant contributions 
followed by several authors [23-26], the latter of which sparked our original interest in the 
subject [27,28]. What particularly attracted our attention was the geometrical setting 
underlying the construction of Polyakov's (world-line) path integral. Within this context 
we pursued the problem of tracing the field theory origins of Polyakov's spin factor, 
introduced by him to properly account for the propagation of a free, spin-1/2 particle- 
like entity on a closed (Euclidean) space-time contour. In Ref. [28] we established, 
via a well-defined procedure, the emergence of the spin factor through the recasting of 
the, spin-1/2, matter field sector of a gauge field theory from a functional to a (world- 
line) path-integral, entering as an appropriate weight to account for spin. At the same 
time, the field theoretical interaction termP] ipj^ipA' 1 is replaced by a 'factorized' Wilson 
line (or loop) which accounts for the effect of the gauge field on the world-line paths, 
equivalently describes its interaction with the matter particles. Our first applications 
turned in the direction of considering situations when it is justified to set the spin factor 
to unity -an occurrence which facilitates a factorization of the infra-red sector of the 
gauge field theory, in its perturbative version [29]. 

More recently, we have tested the possible merits stemming from the aformentioned 
disentanglement between spin-factor and Wilson line(loop), inherent in the Polyakov 

1 The interaction term is displayed generically, independent as to whether the gauge theory is abelian 
or non-abelian. 
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(world-line) path integral for spin-1/2 particle-like entities, as far as the task of facili- 
tating effective, perturbative computations in QCD is concerned [30,31]. In the first of 
these papers the emphasis was placed on extending the world-line methodology to open 
fermionic lines. At the same time we established a procedure by which Strassler's path 
integral expression for spin-1/2 matter particles, entering an (interacting) gauge field 
theory and which contains the term a ■ F in the action, can be recast into the Polyakov 
form which carries, in its place, the spin-factor. Our subsequent manipulations were 
expedited by the presence of the spin factor and produced, as a bonus, the following 
physical picture: Non-trivial spin-factor contributions come precisely at those points 
where a gauge field quantum is emitted or absorbed by the fermionic world-line path. 
Moreover, each such occurence signifies the presence of a derivative discontinuity on the 
path as a four-momentum is locally injected. This is, indeed, a nice intuitive picture 
as it connects the mathematical fact pertaining to the dominance of non-differentiable 
paths on the one hand, with the physical occurence of emission and absorption of quanta 
on the other. In the second paper we focused our attention on the more pragmatic goal 
of developing algorithms, always for spin-1/2 particle (open) world-lines, which lead to 
efficient perturbative computations in QCD pertaining to Green's functions and ampli- 
tudes. Two different alternatives were arrived at according to whether the execution of 
the particle path-integral preceeds or follows the considerations involving the Wilson line 
(loop): (a) The Feynman diagrammatic logic is directly visible and comprehesively dealt 
with, (b) a novel organization of the perturbative expansion is achieved which retains 
the space-time description all the way. Strategy (a) leads to a neat organization of the 
resulting perturbative expression which reduces the computation to two straightforward 
steps. The first pertains, basically, to the spin-factor and amounts to a simple integra- 
tion over a set of Grassmann variables, as many in number as the perturbative order 
considered. The second integration is over a set of proper time parameters inherited 
from the expansion of the Wilson exponential. Possible practical merits of alternative 
(b) constitute an open issue. 

Encouraged by the fact that the Polyakov path- integral for spin-1/2 particle enti- 
ties leads to computational procedures which, both logistically and intuitively, seem to 
present advantages of their own, we undertake, in the present paper, the task of ex- 
tending the application of the relevant methodology to the gluonic sector of QCD. Such 
an effort entails, among other things, the determination of the spin factor pertaining 
to the propagation of a spin-1 particle like entity. This issue will be confronted in Sec- 
tion 2 where we consider the pure gauge field sector of a Yang-Mills system and utilize 
techniques associated with the background gauge field fixing procedure. Focusing on 
effective action terms at the one gluon loop level we proceed, in Section 3, to produce 
a master expression furnishing the Mth perturbative order contribution. The overall 
structure of these terms corresponds to a gluonic world-line loop on which 'vertex oper- 
ators', in the form of plane waves are attached. As in our previous work [31], pertaining 
to open fermionic world-lines, the overall calculation reduces to an integration over a 
Grassmann set of parameters followed by one over a set of Feynman-type parameters. 
The number of variables for each of the two sets is fixed by the perturbative order, while 
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the execution of the integrals themselves is a fairly straightforward matter, as witnessed 
by the applications worked out in Section 4. There, the master formulas are used in 
order to compute the divergent parts of the 2nd, 3rd and 4th order one gluon loop con- 
tributions to the effective action. From the specific manipulations it becomes totally 
apparent how the world-line configuration accommodates, to a given order, the totality 
of the contributing Feynman diagrams via straightforward permutations of the param- 
eters. As already mentioned, our analytic computations in this section are restricted 
to the divergent contributions (furnishing renormalization factors) which can be readily 
identified and isolated from the full expression to each given order. In this connection, it 
can be easily surmised from our master expressions that no divergent terms make their 
appearance above the 4th order, an occurrence that complies with the renormalizability 
of the theory. As far as the finite contributions are concerned, what we can say at this 
stage is that we are in the process of finalizing tests of relevant numerical procedures 
that have been devised for coping directly with their computation. We intend to report 
on this matter in the near future [32]. Finally, in Section 5 we summarize our findings 
and formulate our conclusions, while in an Appendix we trace the main steps involved 
in bringing the spin-factor to its final, ready to apply form. 

2 Polyakov World- Line Path Integral for the Gluon 
Sector of QCD 

The successful trancription of the fermionic sector of a gauge field theory into its Polyakov 
path integral^ form utilizes the fact that the corresponding functional integral is of a 
gaussian type [28]. For the gluon sector, of course, such is not the case. Following Refs. 
[13,19] we proceed by employing the background gauge fixing procedure according to 
which the gauge field splits into a dynamical component, to be denoted by a M , and 
a background field B^. Given that we shall restrict, in the present paper, ourselves to 
the computation of effective action terms, the background field will be considered as 
classical. Let us finally mention that we shall keep our formalism Euclidean throughout 
our analysis. Transcription of our final results to Minkowski space-time will be made in 
the end. In this respect, characterizations such as 'Lorentz generators', 'Lorentz trace', 
etc. will be employed by abuse of language. 

The quadratic part of the (pure) gauge field action reads 

al + c a [(D 2 ) ab ]c\ (2.1) 

where D ab = D ab (B) = d II 5 ab + gf abc B^ is the covariant field derivative in the adjoint 
representation, while c, c are the ghost fields. Obviously, entering Eq. (2.1) is the 

2 For convenience we drop the characterization 'worldline' from hereon, even though we recognize the 
fact that we are using a term that has been established for the characterization of the same author's 
path integral pertaining to string quantization. 
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ab 



{D^D u ) ab - igR 



flU 



4 



Maxwell tensor for the background gauge field, i.e 

K = F;l{B) = -ir bc F^(B) = -if abc (d,Bl - - gr de BlBl). (2.2) 

Adopting the Feynman gauge (£ = f) and introducing the Lorentz generators under 
which four-vectors transform, namely 

( J P°)»v = K S PH S <n> - <WW), (2.3) 

we rewrite Eq. (2.1) as follows 



1 



S 2 = -a; [-(D 2 ) a \ u - g{J ■ F)%] a b u + c a [{D 2 ) ah ] c b 



(2.4) 



In the one loop approximation, to which we shall restrict our considerations in this 
work, the effective action, as a functional of the background field, is given by 

r, |B] = iTVln (-D* — gj ■ F) — TrH-D>) = [B] + [B] . (2.5) 

In what follows it suffices to work with i\ s ; uons [B] as T g h osts [B] is simply given by the 
first of the two terms entering the gluon contribution to the effective action, multiplied 
by (-2). 

Employing Schwinger's parametrization formula [22] we write (trace on 'Lorentz' and 
color indices) 



oo 

1 f dT f 

T^iuons [B] = -- J — J d D xTrK(x,x;T), 



(2.6) 



where 



(2.7) 



corresponds to the (dynamical) gauge field propagator kernel in the background field. 

The world-line path integral for K(y,x\T) ab u results through standard procedures 
(see, e.g., Ref. [30]) and reads 



K(y,x-T)f„ 



I 



x(0)=x,x(T)=y 





T 

~ J dtx\t) 




Dx{t) exp 


Pexp 


y 








-i ab 



ig J dtx-B + g J dtJ-F 
o o 

(2.8) 



-I iiu 



As already established in Refs. [30,31] the Polyakov path integral results once we 
apply the 'area derivative' operator [33,34] given by 



5 



lim dss 



s 2 



Sx^t + %)5x v (t - %) 



(2.9) 
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and use, at the same time, the identities 

T T 
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and 



Pexp(ig J dtx ■ B) ab = Pexp(ig J dtx ■ B)™ 1 {-igF{x{t))f p l a2 P exp(ig J dtx ■ B) 
o * o 

(2.10) 

T T T T 

J dt — exp(-^ J dtx 2 ) = ^ J dtu^ [x(t]] exp(-^ J dtx 2 ), (2.11) 



where uj^ v expresses the rotation of the vector tangent to the trajectory [26] and for 
paths described by differentiable functions assumes the form 



T 



(2.12) 



A more careful discussion pertaining to the spin factor is conducted in the appendix. 

Once performing a partial integration, Eq. (2.8) assumes its Polyakov path-integral 
form which reads 



K(y,x;T)f u 



ii i 
J Dx(t)exp(~ J dtx 2 )Pexp(^ J dtJ ■ lu) ^P exp(ig J dtx ■ B) c 



x(0)=x,x(T)=y 



(2.13) 



In turn, the corresponding expression for the effective action, including the contribution 
from the ghost term, becomes 



ri [B] 



oo 

1 f dT 



T 

x{0)=x(T) 



Dx(t) exp(-i J dtx 2 ) {Tr L $ [1] [x] - 2} Tr c Pexp(ig J dtx ■ B), 

o o 

(2.14) 

where the indices on the traces stand for 'Lorentz' (L) and 'color' (c) and where 



$W [x\ = Pexp 



dtJ-uj[x(t)} 



(2.15) 



-I flU 



is the spin- factor expression^ for the spin-1 particle-like entity (gluon) propagating in 
(Euclidean)space-time. It is not difficult to see that the spin-factor has a restricted 

3 Strictly speaking, Polyakov's spin-factor is given by the (Lorentz) trace of $W 
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dependence on a path's profile. As argued in Ref. [30] and further deliberated on in 
the appendix, contributions of the spin-factor to the path integral comes solely from 
points where a four-momentum is applied through an emission or absorbtion of a gauge 
field quantum. Roughly speaking, this has to do with the fact that the expectation 
value < x fM x u > — < x^x u >, as computed through the path integral, vanishes unless a 
four-momentum is imparted at the point x. 

For the sake of comparison we give the corresponding expression for the one fermionic 
loop expression which contributes to the effective action [28] . It reads (color matrices in 
the fundamental representation) 

oo T T 

V 1Jermi(ms \ B ] = \f d Y f Dx(t) exp(-^ J dtx 2 )Tr L & 1 ^ [x] Tr c Pexp(ig J dtx ■ B) 

x(0)=x(T) 

with the spin factor now given by 



(2.16) 



[x] = Pexp 



dtS-u[x(t)] 



(2.17) 



and where the corresponding Lorentz generators belong to the spinor representation, i.e. 



^ [7^5 iA 



(2.18) 



Generally put, the Polyakov path integral recasting of a relativistic, quantum field theo- 
retical system provides a unified basis for the description of the propagating particle-like 
entity; simply one has to adjust the weight provided by the spin-factor to its particular 
form. Thus, for a spin-zero particle the relevant weight factor is, simply, unity (note, 
in this regard, that ghosts fall into this class irrespective of the anticommutation rela- 
tions, cf. minus sign, they obey) while for spin-1/2 and spin-1 particle-like entities the 
corresponding weights are provided by Eqs. (2.17) and (2.15), respectively. 

For completeness let us mention that the path integral expression for the gluonic 
Green's function, namely 



iG(y,x)f u = I dTK(y,x;T)%„ 



(2.19) 



is determined once the substitution from Eq. (2.13) is made for the propagation kernel. 



3 The one Gluon Loop, M-point Effective Action 

In this Section we shall perform a number of manipulations through which we shall arrive 
to ready-to-apply master expressions for the computation of one-loop effective action 
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terms. Let us commence our efforts by substituting plane wave solutions of the (classical) 
Yang-Mills equations satisfied by the background field, i.e. setting B^(x) = t^ e™e iPn ' x for 
each gauge field entering the M-th order term in the expansion of the Wilson exponential 
in (2.14). We obtain 



oo 

rM(p 1 ,...,p M ) = -±(ig) u Tr c (t?...®) J f 



n=M 



n dt 



9(t M , ti) 



x 



1 x 1 M 

/ Dx(t) Yl e n ■ x(t n ) {Tr L $W [x] - 2} exp -- / dtx 2 + i^p n - x(t n ) 

J rrjn n=M { n=l 



x(0)=x(T) 

+permutations, 



(3.1) 



where 9(tM, ■■■,ti) = [ [ 0(t n +i — t n ) and where the indication permutations refers 

n=M-l 

to all possible rearrangements of the t n and the tf-f associated with the, M in number, 
background gauge fields. 

Our computational strategy for confronting the above quantity coincides with the 
one employed in Ref. [31]. It relies on a move to recast the spin-factor expression into 
an explicitly path-independent form. Once this is done the path integration can be 
immediately performed given that the 'action functional' is a simple gaussian (with a 
linear term). Subsequently, we shall deal with the spin-factor. 

Following the procedure employed in the aformentioned reference we introduce the 
Grassmann variables £ n and £ n through which the e n ■ x(t n ) factors in Eq. (3.1) are 
elevated into exponentials according to [13] 



ie n ■ x(t n ) = J d£ n d£ n exp [i£ n £ n e n ■ x(t n )] 
After substituting in (3.1) we obtain 



(3.2) 



rW(p 1 ,...,p M ) = --g M Tr c (1%'...1%: 



dT 
~T 



n=M 



} I / dtn 



9(tM, h) 



n / di n 



d^n 



X 



Dx(t) {Tr L $ [1] [x\ - 2} cxp 



x(0)=x(T) 

+per "mutations , 



I f M ~ 

-- / dtX 2 + i ^ Htn) ■ x (tn 
n n=l 



(3.3) 



having set 



d 



kfjiitn) — Pn,fj, + CnCn^fxg^ ■ 



(3.4) 
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Recalling the original specification of the spin-factor, which utilizes the employment 
of the area derivative, let us rewrite Eq. (3.3) in the specific form which takes into account 
the fact that the gauge potentials entering the expansion of the Wilson exponential are 
plane waves. For the Mth order term we write 



Trr$ [1] 



M 



exp 



x exp 



^ ] k(t n ) ■ x(t r 



n=l 

M 



Tr^Pexp I —% I dtJ ■ 



5s 



^2k(t n ) ■ x(t n ) 



n=l 



(3.5) 



The above expression illustrates in an immediate, albeit formal, manner the path- 
independence of the spin-factor: One observes that the area derivative acting on the 
exponential will produce delta functions entering the parametric integration, cf. Eq. 
(2.9), entering the definition of the area derivative operator. A well defined argument 
leading to this result is provided in the appendix. 
Returning to the case in hand, we write 



r[ M \p u ..,p M ) = --g M Tr c {t^ ...t%) 



dT 
~T 



n=M 



n dt 



9{tM-, t\) 



n/ 

n=M 



d^ndin 



X 



f 1 } M 

jrr L $ [1] k -2| / Dx(t)exp — dtx 2 + i^fc(i„) • x(t n ) 



x(0)=x(T) 



+per 'mutations . 



(3.6) 



The first task we shall carry out is to perform the, basically Gaussian, path integral. 
Straightforward manipulations, partly displayed in the appendix, lead to the result 



/ 



Dx(t) exp 



x(0)=x(T) 



, y m m 

- J dtx 2 (t)+iYk(t n ).x(t n ) = (27r) D ^(E^)(4^p7^ 



x exp 



' Pm,9nG(t n , t Ti 



n<m 



c-. ^ ^ £,n£,n£,m£,m£ ' £ 9 n d m G(t n , t Tl 



n=£m 



(3.7) 



One notes contributions pertaining solely to the points of attachment of external gauge 
field on the loop contour. In the above expression the following Green's function [13] 
has been employed 



G(t,t') = \t-t'\ 



\t-t' 
T 



(3.8) 
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It corresponds^ to the motion of a one- dimensional particle moving on a closed contour, 
i.e. 



Id 2 1 
-^-G(t,t') = -5(t-t') + - 



(3.9) 



and obeys the boundary conditions G(Q,t') = G(T,t') and G(0, t') = G(T,t'). 

Introducing the dimensionless parameters Ui according to ti = Tv,i,i = l,...,n, the 
interim result for t[ M \ Pi , ...,p M ) reads 



r[ M \ Pl ,..., PM ) = - 



M 00 

1 -9 M (27t) D 6^(f2Pn)Trc(t a G M ■ • ■ ^)^yjf 2 1 dTT 



M-D/2-X 



X 



T ^2 Pn ' PrnG{u n , U r 



Y\ / du n 6(u M ,--,ui)F ( - M) (ui,... } u M ;T)exp 

n=M q L n<m 

+permutati ons , 

where G(u n ,u m ) = \u n — u rn \ [1 — \u n — u m \] satisfies the additional properties 
d n G(u n , u m ) = G(u n , u m ) = sign{u n - u m ) - 2(u n - u m ) = -G(u m , u n ) 

and 

-d n d m G(u n ,u m ) = dlG(u n ,u m ) = G(u n ,u m ) = 2 [S(u n - u m ) - 1] . 
Finally, in Eq. (3.10) we have set 

i 



(3.10) 
(3.11) 
(3.12) 



F^ M \ Ul ,...,u M ;T) 



n=M 



I [ / din&in 



Trr§ [1] 



k 



x exp 



^ ^ ^n^n^ ' Pm^nG (w n , Urn) 



^ ^ ^n^n^rn^m^ ' £ 9 n d m G(u n , Urn) 



n=£m 



(3.13) 



The spin-factor can now be brought into a ready-to-apply form through a series of 
manipulations that are outlined in the appendix. The following result is arrived at 



k 



Pexp 



M 



n=l 



M 



M 71 2 -1 



n=l 



Yl Yl ^2<T 2 («2)0 Pl «T 1 (ni) + ... 

n2=l ni=l 



(3.14) 



4 To be exact, the solution of Eq. (3.9) includes an arbitrary constant which does not appear in Eq. 
(3.8) on the account that momentum conservation has been imposed. 
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with 



</v(n) = 2Un(e; Pn , u - eX) + ^n + i^n + iUn(e; +i el - e n v +1 e;)S(u n+1 - u n ) (3.15) 

and where we have designated that £m+i — £m+i = 0. 

Two observations of practical interest can be made in connection with the above 
expression for the spin-factor. First, it is clear that the number of terms in the expansion 
of the exponential in Eq. (3.14) terminates at M as the saturation point of the Grassman 
variables is by then reached. Second, the delta function containing term in (3.14) implies 
that for a given ordering there is a contribution from coinciding points, u n and u n+ i 
in this case. This occurrence signifies the presence of a 'four-gluon vertex' which is 
automatically included in a given perturbative calculation, along with the (derivative- 
dependent) 'three-gluon vertices' represented by the first term. One thereby concludes 
that the Mth order perturbative contributions to the effective action are classified, via 
the spin-factor, exclusively by the number of the points of gluon (single or pairwise) 
attachements on the closed world-line contour, in all possible permutations. Accordingly, 
the computation of the M-point effective action term will collect all Mth order, 1PI 
Feynman diagrams. 

We mention in passing that fermionic loop contributions to the effective action easily 
follow by referring to Eqs. (2.16)-(2.18). One simply has to make the substitution 
</>(xr(Jpa-)fiv — ¥ S^ v (f)^ v in Eqs. (3.14) and (3.15) and use the fundamental representation 
of the group. 

Given the expressions we have arrived at, what remains to be carried out are the 
integrations over the Grassmann variables as well as the parametric integrations entering 
Eq. (3.10). Numerical methods have been developed for this purpose whose report is 
forthcoming [32]. For the rest of this paper we restrict ourselves to the computation 
of the divergent part of the effective action. In this regard, let us observe, by looking 
at Eq. (3.13), that ultraviolet divergencies will occur only for terms of order M = 
2,3,4. Specifically, by focusing on the terms that have the minimum number of p n ^ 
factors one determines, through dimensional considerations, that they should carry the 
compensating factor T 2 ~ M for M — 2,3, 4. The latter, combines with t m ~ 1 ~ d I' 1 j n (3.10) 
to produce divergent terms ~ T (2 — Further inspection shows that no such terms 
arise for M > 5, a fact that directly complies with the renormalizability of the theory. 

As our final remark in this Section, let us mention that our manipulations up to now 
have treated the gluons as 'real', i.e. we have assumed that 

e n - Pn = and p\ = 0, n=l,2,---. (3.16) 

In the following Section, where our computations will run against infra-red problems, we 
shall protect our expressions from corresponding divergencies by going slightly off mass 
shell, i.e. by setting p 2 n = A 2 , with A 2 > A.q CD . 
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4 Computation of Divergent one-loop Effective Ac- 
tion Terms to Fourth Order 



In this Section we shall apply our comprehesive formulas given by Eqs. (3.10) and (3.13)- 
(3.15) towards the computation of the divergent contributions to the M = 2,3,4 terms 
in the expansion of the effective action -in fact, the only terms which exhibit ultra-violet 
divergencies. We leave the task of computing finite contributions, to the same order, to 
a future paper where numerical methods will be applied. 

4.1 The two gluon contribution (M=2) to the Effective Action 

The present calculation pertains to the one gluon loop configuration with two (trun- 
cated) gluon and ghost attachments. Our master expression accomodates all the con- 
tributing Feynman diagrams. From Eqs. (3.14) and (3.15) we determine, for M — 2, 

Tr L $W = D — ■ e 2 pi ■ p 2 . (4.1) 

Upon substituting in (3.13) and performing the Grassmann integrations we obtain 

F^ M = 2 \ Ul ,u 2 ;T) = ~(D - 2)e l ■ e 2 G( Ul ,u 2 ) - 8s 1 ■ e 2 Pl ■ p 2 . (4.2) 

One notes that the delta function entering the specification of G(u±, u 2 ) accomodates the 
contribution coming from the class of Feynman diagrams wherein the two (truncated) 
external gluons attach themselves to the loop through a four-point vertex. 

The above result when substituted in (3.13) gives, after an integration by parts which 
results in the replacement G(ui,u 2 ) — > — Tpi ■ p 2 G 2 (ui,u 2 ), 

r[ M=2) ( Pl ,p 2 ) = -\{2^\ Pl + P2)Tr c {t a St a S)j^ m - 2 e l ■ e 2 Pl ■ p 2 



x 



OO 1 U2 

J dTT^ D/2 J du 2 J d Ul (D - 2)G 2 ( Ul , u 2 ) - 8] exp [-T\ 2 G{u u u 2 )] 



■^permutations , (4.3) 

where the infra-red cutoff A has been introduced by going off shell. The integrations in 
the last equation can be easily performed and lead to the final result 

i 2 / 2\ 2 -' D / 2 

T[ M=2) (PuP 2 ) = -^r^H^+P^^^U^j e l -e 2 Pl -p 2 

( D\ 11 -7(2- D/2) (D D \ 
XF \ 2 ~ 2) 3 - 2(2 - D/2) B U ' 27 ' { ] 

where the adjustment g 2 — > g 2 D = g 2 fi 4 ~ D was made in order to restore dimensional 
consistency. The term permutations in Eq. (3.10) has been duly taken care of by taking 
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into account all the rearrangements of indices (1,2) and dividing by 2! in order to comply 
with boson non-distinguishability. 

From Eq. (4.4) we verify, once returning to Minkowski space-time, the well known 
result (which does not take into account the contribution from the fermionic loop) 



1 9 2 ,41 
— N— ■ 

2 (4tt) 2 3 2-5/2 



^i-oTT^-Ti— (4-5) 



It is of interest to observe, by referring to Eqs. (2.13), (2.19) and as has been explicitly 
demonstrated in Refs. [30,31], that the corresponding formulas resulting from Polyakov's 
path-integral for open lines have the same basic structure with the ones that have resulted 
from the present considerations pertaining to loops. It, then, becomes a straight forward 
matter to surmise the validity of Ward's identity Z 1 ^ 2 Z g = 1 which is known to hold in 
the framework of the background gauge fixing method. 

4.2 The three gluon contribution (M=3) to the Effective Action 

We now turn our attention to y[ M ~ 3 ^ which summarizes the contributions from the 
classes of Feynman diagrams involving three (truncated) external guons, as well as ghosts. 
Again, our first task is to compute the corresponding expression for the spin-factor. Eqs. 
(3.14) and (3.15) now give 

Tr L $W = D + 8666£3(£ 2 • Pse 3 • p 2 - e 2 ■ e 3 p 2 ■ p 3 ) + 1 • p 3 e 3 ■ Pl - e 1 ■ e 3 Pl ■ j. 

+8£i£i66(^ 1 • P2e 2 ■ pi - e 1 ■ e 2 p 1 ■ p 2 ) 
16 - 

+(e x • e 3 e 2 ■ p 3 - e 2 ■ e 3 e l ■ p 3 )5(u 2 - Ul ) + Te 1 ■ p 3 e 2 ■ Pl e 3 ■ p 2 \. (4. 

The integration over the Grassmann variables can be systematically performed, yield- 
ing the result 

f} m=3 \u 1: u 2 , u 3 ; T) = -^-^{e 1 ■ s 2 [s 3 ■ Pl G(u 3 , Ul ) + s 3 ■ p 2 G(u 3 , u 2 )]G(u l} u 2 ) + 

+S 1 ■ e 3 [e 2 ■ piG(u 2 , ui) + e 2 ■ p 3 G(u 2 , u 3 )]G(u 1 , u 3 ) + e 2 ■ e 3 ^ 1 ■ p 2 G(u u u 2 ) + e 1 ■ p 3 G{ui 
16 

xG(u 2 , u 3 )} + —(e 1 - e 2 e 3 ■ p 1 - e 1 ■ e 3 e 2 ■ Pi)5{u 3 - u 2 ) 
16 

+— (e 1 ■ e 3 e 2 ■ p 3 - e 2 ■ e 3 e l ■ p 3 )5(u 2 - m) + f.t., 

where f.t. stands for 'terms with finite contribution'. Obviously the latter involve terms 
with T to the 0th power or higher, equivalently, they involve more than one (external) 
momentum variables. Let us reiterate that the finite terms should be computable through 
numerical methods that are currently being developed. 
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Substituting the above result in Eq. (3.10) we obtain 

oo 

rS M=3) (pi,p 2 , P3 ) = -i(27r) 4 5( 4 )( Pl +p 2 +p3)Tr c (^t-^0^^ J dTT 1 '^ 

o 

1 U3 U2 

J du 3 J du 2 J dui x {A(D - 2)[e l ■ e 2 e 3 ■ p 2 {u 2 - ui) + e 1 ■ e 3 e 2 ■ pi(l - (u 3 - ui))] 



+e 2 ■ e 3 e 1 ■ p 3 {u 3 - u 2 )\ - l^e 1 ■ e 2 e 3 ■ p 2 + e 1 ■ e 3 e 2 ■ pi)5{u 3 - u 2 ) 
-16(e 1 .e 3 e 2 ■ Pi + e 2 ■ e 3 e 1 ■ p 3 )5(u 2 - Ul ) + f.t.} 
TX 2 

X exp{ ^-[( M 2 - «l)(l - (U 2 ~ Ml)) + («3 - M 2)(l - («3 ~ "2)) 

+ (u 3 — ui)(l — (u 3 — Mi))]} + permutations. (4.8) 

It is easy to see that the first term in the curly brackets takes care of the Feynman dia- 
grams involving three-gluon vertices while the other two, which carry the delta-functions, 
collect the contributions from diagrams with one four-vertex. To further guide the reader 
let us also mention that use was made of Eqs. (3.11) and (3.12). Accordingly, the above 
expression refers to the specific ordering which enters these equations and underlies the 
particular integrations over the parameters u±, u 2 and u 3 , an occurrence that will be 
rectified shortly. Finally, in the exponential factor we have set 

v\=vl=v\ = A 2 , 2pi • p 2 = 2 Pl .p 3 = 2p 2 ■ p 3 = -A 2 . (4.9) 

Next, we make the variable change u 2 — u\ = x 2 and u 3 — u\ = x 3 which casts Eq. (4.8) 
into the form 

1 2 / 2\2-D/2 

r[ M = 3 \ Pl ,p 2 ,p 3 ) = --(27r) 4 ^)( Pl + P2 +Pz)Tr c {t a St^) 9j ^ 

x{-A(D - 2)a D (e 1 ■ e 2 e 3 ■ p 2 + e 1 ■ e 3 e 2 -pt + e 2 - e 3 e 1 ■ p 3 ) + 8b D (e l ■ e 2 e 3 ■ p 2 

+3e x ■ e 3 e 2 ■ pi + 2e 2 ■ e 3 e 1 ■ p 3 )}T(2 - — ) + f.t + permutations, (4.10) 

2 

where we have introduced 

1 x 3 

a D = 2 2 ~ D/2 J dx 3 J dx 2 x 2 [x 2 (l - x 2 ) + (x 3 - x 2 )(l - x 3 + x 2 ) + x 3 (l - x 3 )] D/2 ~ 2 
o o 

(4.11) 

and 

i 

b D = J dx 3 (x 3 (l-x 3 )) D/2 ' 2 . (4.12) 
o 
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Obviously — £ and 64 = 1. 

In order to obtain the final result we need to take into account contributions coming 
from all permutations of the variables U\, u 2 , u 3 and divide by 3! to compensate for 
boson indistinguishability. The result can be easily obtained using the cyclic properties 
of the trace. One, finally, obtains 

1 2 / 2\2-D/2 

rr =3 V,P2,P3) = -(2n) 4 S^( Pl + p 2 +Ps)Tr c (t a Jt a G H^gj^ ^tt^J 
x{A{D - 2)a D {e 1 ■ e 2 e 3 ■ p 2 + e l ■ e 3 e 2 ■ pi + e 2 ■ e 3 e 1 ■ p 3 ) - m D (e l ■ e 2 e 3 ■ p 2 
+6 1 ■ e 3 e 2 ■ Pl + e 2 ■ e 3 e 1 ■ p 3 )}T(2 - y ) + /.*. (4.13) 
in agreement with the known result. 

4.3 The four gluon contribution (M=4) to the Effective Action 

The computation in the present subsection pertains to a 'world-line' diagram which 
collectively accomodates all the classes of the contributing Feynman diagrams, i.e. with 
zero, one and two four- vertices (plus, of course, contributions from ghost diagrams). As 
our present analytic computations refer to the divergent part, let us isolate the relevant 
contribution (terms with the factor 1/T 2 ) entering the expression for the spin-factor. We 
find 

32 - 

Tr L $W — D + —UUs^teMe 1 ■ e 4 e 2 ■ e 3 - e 1 ■ e 3 e 2 ■ e 4 )5{u, - u 3 )5(u 2 - Wl ) + f.t. 

(4.14) 

Integration over the Grassmann variables is a straightforward matter and gives 
F iM =^ ( U4 , u 3 , U2, u 1 ;T) = [e 1 ■ e 2 e 3 ■ e 4 G(u u u 2 )G(u 3 , u 4 ) 

+S 1 ■ e 3 e 2 ■ e 4 G(ui, u 3 )G(u 2 , u 4 ) + +E 1 ■ e 4 e 2 ■ e 3 )G(ui, u 4 )G(u 2 , u 3 )] 
+^{e 1 ■ e A e 2 ■ e 3 - e l ■ e 3 e 2 ■ e 4 )5(u 4 - u 3 )5(u 2 - Ul ) + f.t. (4.15) 

Substituting the above expression into Eq. (3.10) we get 

r[ M=4 \puP2,P3,P,) = ~(2n) 4 S^( Pl +p 2 + p 3 + p 4 )Tr c (t a G H a G n a Jt a G ^ 9 ° 



x 



2^ " 1 ^ z 1 < fi>^ ■ ov-G "G "G "G J (27r) D / 2 

1 U4 M3 U2 

J du 4 J du 3 J du 2 J d Ul {A(D - 2)(e 1 ■ e 2 e 3 ■ e 4 + e 1 ■ e 3 e 2 ■ e 4 + e 1 ■ e 4 e 2 ■ e 3 ) 
0000 

.1 Jl?> „4rf/„, „, \ . p/ , \1 Al-r, r>\ 1 4 2 3 s 



-4(D - 2)e 1 ■ e 2 e 3 ■ e 4 [5(u 2 - Ul ) + 5(u 4 - u 3 )\ - 4(D - 2)e 1 ■ e 4 e 2 ■ e 3 5(u 3 - u 2 ) + 
+4(D - 2)e 1 ■ e 2 e 3 ■ e 4 5(u 2 - Ul )5(u 4 - u 3 ) + 32(e 1 • e 4 e 2 ■ e 3 - e 1 ■ e 3 e 2 ■ e 4 ) 

4 4 

^2 Pn ' PmG(u n ,U m ) 

_n=l m=n+l 

+ f.t. + permutations. (4.16) 



x5(u 2 - ui)5(u A - u 3 )} 



D/2-2 
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One can easily verify that the first term inside the curly brackets represents contributions 
corresponding to the Feynman diagrams with no four-vertices, the next two to those 
with one and the last to those with two. Of course, the above expression pertains 
to a particular ordering of the variables ui, u^, u%, U4 as reflected in the explicit delta 
functions which make their entrance. 

Performing the parametric integrations, in the specific ordering that appears in Eq. 
(4.16), one obtains 

r[ M=4 \p 1 , P2 ,Ps,p 4 ) = ~(2n) 4 S^( Pl +p 2 +p 3 +Pt)Tr c {1%1%1%t%)-J^ 

x{A(D - 2)A D {e l ■ e 2 e 3 ■ e 4 + e 1 ■ e 3 e 2 ■ e 4 + e l ■ e 4 e 2 ■ e 3 ) - 4(D - 2){B D - C D )e 1 ■ e 2 e 3 ■ e 4 - 



D 



-A(D - 2) J D D e 1 • e 4 e 2 ■ e 3 + 32C D (e l ■ e 4 e 2 • e 6 - e 1 ■ e 3 e 2 ■ e 4 )}T(2 - —) + f.t. 
+permutations, 

where we have set 



(4.17) 



A 



D 



1 U4 U3 U2 

dui J du 3 J du2 J du\ 
oooo 

1 U4 U3 «2 



4 4 



^2 Y Pn ' PmG(u n , U m ) 



,n=l m=n+l 



D/2-2 



B 



D 



c D = 



J du4 J ditz J du 2 J dui[5(ii2 — Ui) + 5(ii4 — Uz)] 



U4 U3 U2 

du4 / dus / du2 / dui5(u2 — Ui)5(ii4 — 113) 



4 4 



^2 Y Pn' P m G(u n ,U m ) 
n=l m=n+l 



D D = 



0000 

1 «4 U3 U2 

du4 J du 3 J du 2 J dui5(u 3 — u 2 ) 
0000 



4 4 



4 4 



Y Y2 Pn ' VmG{u n , U r 
D/2-2 



,n=l m=n+l 



D/2-2 



id U m ) 



n=l m=n+l 



(4.18) 



One trivially finds A4 — |, B4 — |, C4 = \ and D4 — |. 

The remaing step is to perform all reorderings of the u variables and divide by 4!. In 
this way one arrives at the final expression 

r[ M=A \pi,P2,P3,P,) = -\(2n) 4 S^( Pl +p 2 +p, + Pi )Tr c (t a G %H a ^)j^^ 

x {4(D - 2)A D (e 1 ■ e 2 e 3 ■ e 4 + e 1 ■ e 3 e 2 ■ e 4 + e 1 ■ e 4 e 2 ■ e 3 ) - 2(D - 2) 

x{B D -C D + Do^e 1 ■ e 2 e 3 ■ e 4 + e 1 ■ e 4 e 2 ■ e 3 ) + WC^e 1 ■ e 2 e 3 ■ e 4 + e 1 ■ e 4 e 2 ■ e 3 ) 



32C D e 1 -e 3 e 2 -e 4 }T(2-^) + f.t., 



(4.19) 



which is in full agreement with the known results. 
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5 Concluding comments 



Given the schemes pioneered by Bern and Kosower [6] and reformulated by Strassler [13] 
based on string and world-line agents, respectively, and which aim at expediting pertur- 
bative computations in QCD both economically and efficiently, it becomes important to 
assess the relevant merits of yet another competitive proposal advanced in the present 
paper which utilizes the Polyakov world-line path integral. Directing, to begin with, our 
comments towards making comparisons with Strassler's approach we could say that the 
basic difference between the two world-line based schemes is how the disentanglement, 
between the weight factor pertaining to the spin of the propagating particle-like object 
on a given path and the dynamical factor represented by the Wilson line (loop), is ac- 
complished. In Strassler's case this task is confronted by using super-particle degrees 
of freedom (one dimensional) and generates a term in the corresponding Lagrangian of 
the form ■ip fl F tlu if) u . In the Polyakov (world-line) version, on the other hand, the issue 
is addressed via the introduction of the spin-factor. We believe that the separation, 
featured by the latter scheme, between 'geometrical' characteristics of paths on the one 
hand and dynamics -as embodied in the Wilson line (loop) factor- on the other, leads to 
an organization of the path integral expression which further facilitates the 'efficiency 
factor' for performing perturbative computations. In particular, it offers a unified base 
for treating, spinors, gauge fields and ghosts; all one has to do is adjust the master for- 
mula, which yields the computational rules, to the appropriate spin factor. Moreover, 
it lends itself to straightforward extensions for applications to processes involving open 
fermionic world-lines, as established in Refs [30, 31]. Referring, finally, to the string- 
based approach of Bern and Kosower we note that the pinching issue, which arises in 
the non-abelian case and whose confrontation requires the application of a certain set of 
mnemonic rules, does not arise in our approach. We expect to further demostrate the 
virtues of the Polyakov world-line path integral scheme toward the calculation of two 
gluon loop contributions to the effective action by generating the corresponding master 
formulas. 
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Appendix 



In this Appendix we shall pay closer attention to the spin factor with respect to both 
carrying out the path integral in Eq. (3.3) and establishing the result encoded in Eqs. 
(3.14) and (3.15). Looking at the identity given by Eq. (2.11) we present the proper 
(regularized) expression for the tensor reads as follows 



T T 



J dtajn„[x(t)] = lim^ J ds J dt 2 J dt^x^^x^) - Xv^x^t^Sfo - h - s). 

(A.l) 



-e 



Now, if the functions (on the line) x^{t) are infinitely differentiable, then we can, once 
taking into account that \t 2 — 1±\ < e, write i M (t 2 ) = x^ti) + 0(s) as well as x u {t 1 ) = 
Xu(h) + CK S ) an d immediately conclude that 



r 



J dtu^[x(t)} = ^ dt[x^t)x u (t) -x u (t)x^(t)}. (A.2) 



Otherwise, one should use the limiting expression according to (A.l) when performing 
manipulations that involve the spin factor. 

Let us proceed with the computation of the path integral entering Eq. (3.3). We set 

l, u = Jd 4 a J Dx(t)^[x(t)], u exp{-S[x]}, (A.3) 



x(0)=x(T)=a 

where 

M 



S[x] = - j dtx 2 (t)-tjrktn)-x(t n ). (A.4) 

n=1 

To compute 1^ we make the variable change x — > x + x d , where x cl is a solution of 
the classical equation of motion resulting from the above action. Specifically, we have 



M M 



x%(t) = -2i ^ k(tn)5(t - t n ) => x*(t) =2iJ2 HQ A(t, Q + a, (A.5) 

n=l n=l 

where we have employed the Green's function 

A(i, if) = t ^Ae{t' - t) + t -^ F ^6(t - t% A(0, t') = A(T, t') = 0. (A.6) 
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The new action functional is now specified by 

T 

S[x] -> ^ J dtx 2 {t) + S[x% (A.7) 
o 

where 

MM M 

S \ xd \ =^2^2kQ ■ Ht m )A(t n ,t m ) -iJ^Pn-a. (A.8) 

n=l m=l n=l 

We immediately observe that integration over a (translational zero modes) leads to 
momentum conservation which enters Eq. (3.7). The rest of the expression for S[x d ] 
produces the terms entering the exponential factor in the same equation. 

Turning our attention to the spin factor we first note that the variable change x — > 
x + x cl leads to 

T T T 

J dtuj^[x}^ J dtco^[x cl ] + ^ J dt[x„(t)x u {t) - x u (t)x^{t)}, (A.9) 



having taken into account that the contours x(t) are to be integrated with respect to a 
quadradtic action functional (cf. Eq. (A.7)), which implies, cf. Eqs. (A.l) and (A. 5), 
that mixed terms in x and x cl drop out. Let us finally note that for paths that are 
infinitely differentiable, in which case Eq. (A. 2) strictly holds true, the integration of 
the spin factor with respect to the quadratic action functional yields unity. All this leads 
to the following result as far as performing the path integral in Eq. (A. 3) is concerned. 

M 



i„ = (j>) j^sj-^H 



x exp 



' Pm,9nGr(t n , t m ) 

n<m n^m 



„ ^ ^ ^n^n^m^mE- ' £ d n m G(t n , t m ) 



2 



(A.10) 



The above result explicitly demonstrates our assertion that the overall contribution from 
the spin factor is exclusively determined by those points on a given path where a mo- 
mentum is imparted via the action of an external gauge field. 

The final result is obtained once we substitute (A. 5) into (A.l). We get 

T . M MM 



n=l n=0 m=0 

e 



d 

x / ds—5(t n -t m -s). (A.ll) 
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The correct handling of the last term follows once we take into consideration that, first, 
m 7^ n on account of the Grassmann variables and, second, it is to be integrated in 
consistence with the time ordering implicit in Eq. (3.1) of the text. Specifically, we have 

T T T e 

...J dt n+ i J dt n J <ft n _ 1 0(i n+1 - f n )0(f n - t n _i) J ds—S(t n -t m -s)...= 





T T T 



= ... J dt n+ i J dt n J dt n -i9(t n+ i — t n )9(t n — t n -i) 



X [28n+l,mfi(tn+l ~ t n ) ~ 2(5^, >rj ,_ 1 (5(t rl — t n -i)]..., 

which allows us to return to Eq. (A. 11) and infer that 



(A.12) 



dtuj^[x cl ] 



"=° 
M 

-^E^^W - e n u +1 e;)5(t n+1 - t n ) (A.13) 



n=l 



With the above results in place, Eqs. (3.14) and (3.15) in the text follow directly. 

The careful course of reasoning we have followed in this appendix can be circumvented 
by the more formal line of procedure adopted in the text. Thus, the validity of the 
aformentioned equations can be established once we observe that 



cxp 



M 



-i J^fc(t n ) -x(t n ) 



n=l 
T 



dt 



exp 



M 



i ^ Htn) ■ X(t r , 



n=l 



M M 



= - Jim / dt^2^2 / dssk fl (t n )5(t n - t - ^)k u (t m )5(t m - t + |) 

n=1 m=1 -s 

T £ 

= - lim I dt I dssxf(t + -)x c J(t - -) 



-e 
e T T 



= lim^ J dt 2 J dt^x^x^h) - x^x^t^Sih-h- s). (A.14) 
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